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The  exponent  0.7  fell  out  from  the  Knuds en.  number  (Kn)  in  an  expres¬ 
sion  occurring  in  equation  39  (page  40)  and  in  the  lower  center  of 
figure  17  (page  41).  Please  add  the  exponent  in  these  two  places. 
The  expression  will  then  correctly  read: 
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ABSTRACT 


Under  contract  to  the  Army  Research  Office,  United  Technology  Center 
has .conducted  an  experimental  investigation  to  determine  the  drag  coefficient 
of  particles  in  flow  regimes  encountered  in  a  rocket  nozzle.  The  acquisition 
of  theoe  data  leads  to  more  reliable  predictions  of  nozzle  performance  ineffi¬ 
ciencies  owing  to  gas -particle  flow. 

''  °o  u 

The  Mach  number- Reynolds  number  regime  traversed  by  a  particle  in  a 
rocket  nozzle  is  described.  The  experiment  to  determine  drag  coefficient 
data  in  this  flow  regime  consists  of  the  electrostatic  acceleration  of  microh- 
sfze  particles  to  sonic  velocities  and  detection  of  their  velocity  decay  in  a 
chamber,  conditioned  to  provide  the  desired  flow  paramete;  s.  The  operation 
of  the  experiment,  method  of  data  reduction,  and  analysis  of  the  experimental 
error  are  presented* 

0  C' 

The  data  are  reduced  in  terms  of  a  nondimens ional  dr?  ,,  coefficient  and 
are  correlated  with  comparable  data  obtained  in  other  flov  regimes.  An 
empirical  relation  is  generated  for  the  drag  coefficient  as  „  function  of 
Reynolds  and  Knudsen  numbers.  This  relation  is  recomm,.;  led  for  use  in 
calculations  of  gas- particle  flow  in  rocket  nozzles. 
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speed  of  sound 
representative  area 
drag  coefficient 

nondimensional  drag  coefficient 

drag  coefficient  for  incompressible  fiow 

"inviscid"  drag  coefficient 

free-molecule-flow  drag  coefficient 

distance  between  Faraday  cases 
total  range  distance 
signal  voltage 

drag  force  i 

current 

Knudsen  number 

number  of  measurement  stations  in  tracking  section 

number  of  measurement  stations  in  velocity-measuring  section 

particle  mass 

Mach  number 

pressure 

statistical  factor 

particle  charge 

particle  radius 

resistance 

Reynolds  number 

speed  ratio 

distance  traveled  by  particle 
gas  temperature 
particle  temperature 
time 

zero  reference  time 

particle  velocity  in  velocity-measuring  section 

gas  velocity 

initial  particle  velocity 

particle  velocity 

accelerating  potential 


*'  Greek 

constant  related  to  drag  coefficient 
ratio  of  specific  heats 
variable 
gas  viscosity 
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o  variance 

Pg  gas  density 

p  particle  density 

\  mean  free  path 
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1.0  INTRODUCTION 


The  use  of  a  metallic  fuel,  i  uch  as  aluminum,  in  rocket  propellants  aug¬ 
ments  the  available  energy  per  u.  It  mass  and  improves  performance.  The 
improvement  is  not  as  large  as  thermodynamic  calculations  indicate  because 
the  micron,  size  particles  of  metal  oxide  produced  upon  combustion  of  the 
metal  and  carried  out  by.  the  exhaust  gases  are  unable  to  maintain  kinetic  and 
thermal  equilibrium  with  the  exhaust  gases ;  this  results  in  a  performance 
inefficiency..  Analyses  of  gas  rpar.ticle  flow  in  rocket  nozzles  have  been 
developed  to  predict  this  inefficiency^*  but  fundamental  data  for  the  aero¬ 
dynamic  drag  force  on  the  particles  have  not  been  available.  The  aim  of  this 
project  was  the  experimental  determination  of  this  aerodynamic  force  in  flow 
regimes  encountered  in  a  rocket  nozzle. 

The  nondimensional  parameter  which  characterizes  the  aerodynamic 


drag  force  on  a  particle  in  .a  gas  str  eam  is  the  drag  coefficient, 
as 


CD 


D 


1  ,  .2  A 
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where  is  the  drag  force,  p  the  gas  density,  -u  the  gas  velocity,  u  the 

particle  velocity,  and  A  a  repre  sentative  particle  area.  In  general,  the  drag 
coefficient  of  an  object  is  dependent  on  several  parameters  such  as  shape  and 
orientation  with  respect  to  the  flow,  M,  Re,  and  turbulence  level.  Experi¬ 
ments^  have  shown  that  particles  found  in  rocket  nozzles  are  spherical 
in  shape,  and  order -of ^magnitude  analyses  indicate  that  the  particle  drag 
coefficient  is  primarily  a  function  of  M  and  Re. 
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The  M-Re  regime  traversed  by  a  particle  in  a  typical  rocket  nozzle  is 
illustrated  in  figure  1.  A  complete  gamut  of  flow  regimes  from  continuum  to 
fre  e -molecule  flow  are  encounter  ed  as  the  particle  moves  from  the  ro.cket 
.chamber,  through  the  throat,  and  into  the  expansion  section.  No  data  .exist for 
the  drag  coefficient  of  a  sphere  iu  the  flow  regimes  bounded  by  Re  less  than 
l.,000  and  M  less  than  2.  A  reasonable  formula  for  the  drag  coefficient  inthis 


*  Parenthetical  superscript  numbers  denote  references  appearing  .on  page  46. 
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region  has  been  dovised(2)  but  its  validity  is  questionable  until  confirmed  by 
experimental  data. 

c 

The  following  sections  describe  the  development  and  operation  of  an  ex¬ 
periment  to  obtain  the  desired  drag  coefficient  data.  The  results  are  dis¬ 
cussed  and  an  empirical  equation  representing  the  data  is  presented. 
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2.  0  DESIGN  AND  DEVELOPMENT  OF 
THE  EXPERIMENTAL  TECHNIQUE 


There  are  two  primary  methods  by  which  the  drag  coefficient  of  an  object 
can  be  measured;  one  is  to  suspend  the  object  in  a  gas  stream  (e.  g. ,  in  a  wind 
tunnel)  and  measure  the  drag  force  directly,  and  the  other  is  to  fire  the  object 
in  a  ballistic  range  and  deduce  the  drag  coefficient  from  time -distance  data.  The 
major  difficulty  in  using  the  wind-tunnel  technique  to  obtain  data  in  the  flow  regime 
of  interest  lies  in  the  requirement  for  accurate  measurement  of  very  small  forces. 
The  design  complexities  of  a  large,,  low -density  wind  tunnel  and  suspension 
system*  together  with  extremely  sensitive  force -measuring  instrumentation, 
-suggest  the  impracticality  of  the  wind-tunnel  technique  for  this  study. 


The  small  drag  forces  that  are  characteristic  ofthe  M-Re  regime  of  in¬ 
terest  also  preclude  the  use  of  a  conventional  ballistic  range.  The  difficulty 
lies  in  the  prohibitively  long  distances  that  a  reasonably  sized  projectile 
would  travel  before  a  detectable  change  in  velocity  occurs.  Employing  Stokes' 
drag  law  to  provide  Indicative  results  and  suggest  the  significant  parameters, 
yields 
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where  S  is  the  distance  traveled  by  a  spherical  projectile  with  radius  r^ 
before' its  velocity  has  decayed  to  u^.  The  projectile's  initial  velocity  is  uQ, 

v' 

and  p_  the  projectile's  density,  and  p  the  gas  viscosity.  The  equation  indi- 

A  '  ”  c  <,  ■  8 

cates  that  a  1  -mm  projectile  with  a  density  of  3  gm/cc  and  fired  in  air  at 
3Q0  m/sec, will  travel  250  m  before  a  10%  velocity  decay  ensues.  The  above 
equation  suggests,  however,  that  the  ballistic  range  method  is  effective  if 
very  small  projectiles  can  be  accelerated  to  a  sufficiently  large  velocity  and 
if  their  time  histories  during  deceleration  can  be  recorded  accurately.  The 
calculations  show  that  the  projectiles  must  be  between  i  micron  and  50  microns 
in  diameter  and  must  achieve  velocities  up  to  300  m/sec  to  provide  data  in 
the  flow  regime  of  interest. 


Acceleration  of  micron? sized  particles  to  high  velocities  has  been  accom¬ 
plished,^  experimental  investigations^*^  of  meteoric  impact.  Briefly,  the 
technique  consisted  of  charging  particles  by  subjecting  them  to  a  large 
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potential  gradient  while  in  contact  witn  a  conducting  surface,  and  then  accel¬ 
erating  these  charged  particles  through  a  large  potential  difference.  Assum¬ 
ing  charge-to-mass  ratios  typical  of  these  experiments,  calculations  show 
that  a  potential  difference  of  100  kv  is  sufficient  for  drag -coefficient  experi¬ 
ments.  The  size  of  the  particles  impacting  on  a  metal  sample  was  determined 
by  measuring  the  charge  and  velocity  of  the  accelerated  particles  and  using 
the  work-energy  equality 

qv  =  \  rnttf2  (3)  ’ 

where  q  is  the  particle  charge,  V  the  accelerating  potential,  m  the  particle’s 
mass,  and  u^the  final  particle  velocity.  The  velocity  and  charge  were  meas¬ 
ured  using  Faraday  cages  which  are  small  tubes  through  which  the  particles 
pass  and  on  which  they  induce  their  charge.  The  Faraday  cages  were  con¬ 
nected  to  the  grid  of  a  cathode-follower  circuit  and  the  particle  passage  was 
detected  as  a  voltage  signal  on  ah  oscilloscope*  This  scheme  was  directly 
adaptable  to  measuring  particle  size  in  the  present  study. 

A  schematic  diagram  of  the  experimental  apparatus  is  given  in  figure  2. 
The  acceleration  section,  differential  pumping  sections,,  and  tracking  section 
were  connected  by  400-micron  orifices  through  which  the  particle  passed  as 
it  proceeded  from  the  acceleration  to  the  tracking  section.  The  acceleration 
section  was  connected  to  a  diffusion  pump,  while  mechanical  pumps  were 
sufficient  for  the  differential  pumping  sections. 

In  order  to  produce  the  range  of  Re  desired,  the  pressure  in  the  tracking 
section  ranged  between  1  and  100  torr.  The  acceleration  section,  however^ 
was  maintained  at  10 -4  torr  for  maximum  electrical  insulation. (5)  Feasibility  ’ 
calculations  indicated  that  this  pressure  differential  could  be  maintained  by 
using  two  differential  pumping  chambers  and  orifices  smaller  than  500  microns 
in  diameter  between  the  four  chambers.  The  differential  pumping' chamber 
adjacent  to  the  acceleration  chamber  served  as  the  space  where  the  charge 
and  velocity  of  the  particle  were  measured  in  order  to  determine  its  size; 

It  was  initially  postulated  that  a  Faraday  cage  tracking  technique  would 
be  ineffective  in  the  tracking  section  because  the  particles  would  rapidly  lose 
their  charge  in  the  higher  pressure  environment.  Thus,  it  was  initially 
decided  to  use  a  tracking  scheme  in  which  the  particles  traversing  a  laser 
beam  would  reflect  light  onto  a  photomultiplier  tube.  However,  a  pilot  experi¬ 
ment  performed  early  in  the  program  demonstrated  that  a  charged  particle  at 
atmospheric  conditions  v/ould  not  lose  its  charge  fast  enough  to  make  the 
Faraday  cage  technique  ineffective.  It  was  also  decided  to  use  a  series  of 


Faraday  cages  to  measure  the  particle’s  time  history  in  the  tracking  section, 
m  lieu  of  the  more  complex  optical  technique.  The  details  of  each  section  o% 
the  apparatus  are  described  below, 

**  *  c 

o 

2. 1  ACCELERATION  CHAMBER  .  .  .  , 

A  diagram  indicating  the  essential  components  of  the  acceleration  cham¬ 
ber  is  shown  in  figure  3.  A  magnetized  needle  was  mounted  on  a  micrometer 
drive  directly  above  a  si  tall  hole  (~80  microns  in  diameter)  in  a  thin  dia¬ 
phragm.  Magnetic  part.cles  aligned  themselves  in  strings  along  the  magnetic 
lines  of  force  frem  the  needle.  A  potential  difference  of  from  15  to  2d  l;v  was 
maintained  between  the  needle  and  the  diaphragm.  As  the  needle  was  advanced 
toward  the  hole  b>  the  micrometer  drive,  the  electrostatic  forces  overcame 
the  magnetic  ones,  and  one  charged  particle  was  removed  from  the  extremity 
of  the  particle  string.  The  charged  particle  proceeded  through  the  hole  and,, 
into  the  lower  part  of  the  chamber  where  it  was  accelerated  through  an  addi-  , 
tional  30  to  60  kv.  A  sharp  conical  section  was  used  for  the  ground  electrode 
to  aid  in  alignment  of  the  particle  trajectory  by  causing  the  electrical  lines 
of  force  to  converge  toward  the  cone.  It  was  found  that  the  sharp- edged  cone 
was  not  sufficient  to  align  the  particle  trajectory;  probably  local  charge  ac¬ 
cumulations  on  the  .plexiglass  walls  of  the  container  disturbed  the  electric 
field  symmetry  and  caused  trajectory  misalignment.  The  problem  was  over¬ 
come  by  locating  a  copper  cylinder  in  the  chamber  coaxial  with  the'  centerline 
of  the  system,  which  ensured  a  uniform  charge  distribution. 

Carbonyl-iron  particles  forming  strings  on  the  tip  of  a  magnetized  needle 
are  shown  in  figure  4.  Generally,  only  one  particle  at  a  time  was  removed 
by  electrostatic  forces.  When  the  supply  of  particles  on  the  needle  became 
depleted,  the  needle  was  withdrawn  into  the<  particle  reservoir  and  acquired  a 
new  supply. 

The  acceleration  chamber  mounted  in  place  can  be  seen  in  figure  5.  It 
was  fabricated  of  plexiglass  to  ensure  sufficient  electrical  insulation  of  the 
high  voltages.  The  copper  cylinder  inserted  to  make  the  electric  field  sym¬ 
metric  is  visible.  The  chamber  was  connected  to  an  oil  diffusion  pump 
through  an  angle  valve  and  liquid-nitrogen  baffle.  Design  calculations  using 
the  manufacturer's  specifications  for  the  speed  of  the  diffusion  purr  p  and  con¬ 
ductance  of  the  baffle  predicted  that  a  pressure  of  10“4  torr  could  be  main¬ 
tained  in  the  acceleration  chamber.  However,  under  actual  operation  no 
pressure  lower  than  1Q"^  torr  was  achievable.  At  this  pressure  nc  serious 
electrical  discharge  problems  were  encountered  at  voltages  up  to  70  kv. 

The  acceleration  chamber  was  designed  for  easy  removal  of  the  top  to 
facilitate  cleaning.  The  removed  top  is  shown  in  figure  6,  The  plate  which 
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Schematic  Diagram  of  Acceleration  Section 
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Figure  4*  Electron-Microscope  Photograph  of  Carbonyl-Iron 
Particles  on  the  Tip  of  a  Magnetized  Needle  , 
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,Sjsppb:?ted  the  thin  diaphragm  with  the  hole  beneath  the  needle  was  suspended 
from  the  plexiglass  top  by  three  metal  bolts,  one  of  which  projected  above 
#io  r.uriaes  to  provide  an  electrical  contact  to  the  plate  and  diaphragm.  The 
that  houses  the  micrometer  feed  and  needle  could  alsv  he  removed  and 
disassembled  L~s  cleaning  and  to  replenish  the  particle  supply.  A  close  fit 
d  tho  tep  ia  tho  acceleration  chamber  was  necessary  to  maintain  alignment  of 
'  tho  icccdlo  v/ith  the  system.  '  . 


A  power  supply  capable  of  delivering  150  watts  at  100  kv  provided  the 
aseoTjST-t'lnj  potential*  The  voltage  divider  Los  which  supplied  the  desired 
yoltegoa  .jmd  served  as  a  load  for  the  power  supply  Is  shown  in  figure  7. 

Gcveti  r  esistor  a  constituted  the  voltage  dividing  network;  four  iOM  resistors 
at  the  top,  two  $0b»  resistors  below,  and  a  precision  i  OOK  resistor  at  the 
fcett&m.  The  lead-in  from  the  power  supply,  at  the  top  right  of  figure  7  was 
connected  to  the  needle  housing  through  a  protective  5QQM  resistor.  The  lead 
f&Dia  the  holt  &•  Ending  above  the  top  of  the  acceleration  chamber  was  con¬ 
nected  through  a  1,P0GM  protective  resistor  to  the  juncture  of  the  10M  resis¬ 
tors  which  provifbSi&ths  desired  potential  difference  between  the  needle  and 
the  .thin  4daphrag$m  The  voltage  across  the  10GK  i  esistor  was  proportional 
to  the  total  accelerating,  potential  and  was  measured  by  the  electrometer 
shown  in  the  lower  righthand  corner  of  figure  5.  Each  component  resistor 
waa.  calibrated  at  operating  potentials  so  that  the  voltage  division  latio  and  the 
accelerating  potential  were  known  to  within  ±\%. 


The  pressure  in  the  acceleration  chamber  was  measured  by  a  Pirani 
gauge  located  in  the  angle  valye  adjacent  to  the  chamber. 


2*2  DIFFERENTIAL,  PUMPING  AND  TRACKING  CHAMBERS 


Two  differential  pumping  chambers  were  used  between  the  acceleration 
and  tracking  chambers.  The  velocity-measuring  chamber,  in  which  particle 
velocity  and  charge  were  determined,  was  adjacent  to  the  acceleration  cham¬ 
ber.  The  chamber  adjacent  to  the  tracking  section  served  only  as  a  differ  - 
cntial -pumping  chamber.  The  four  chambers  were  connected  by  three  400- 
fnicron  Orifices. 


The  velocity  measuring  section,  which  was  fabricated  of  aluminum,  can 
bo  seen  in'  figure  5  directly  beneath., the  acceleration  chamber.  It  is  connected 
to  a  7-Hter/oec  Welch  pomp  through  a  2-in.  -diameter  flexible  hose.  The 
presence  was  sufficiently  low  (5  by  10-3  torr)  so  that  the  particle  encountered 
negligible  aerodynamic  drag  while  traversing  this  section.  Therefore,  *he 
particles-  kinetic  energy  in  this  section  was  exactly  equal  to  the  electrical 
work  done -in  accelerating  the  particle.  The  relation  expressed  by  equation  3 
iao  employed  to  determine  the  mass  of  the  particle. 
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Xiio  v^locity>measuring  section  w«»s  closed  on  one  side  by  a  brass  plate 
wuicli  extended  the  entire  length  of  the  differ  antial  pumping  and  tracking 
Cactsoao,  The  Faraday  cage  system  and  cathode-follower  amplifiers  were 
stated  on  the  brass  plate  which  is  shown  in  figure  5.  The  opposite  side  of 
f'Ji'rrXo  >: over ed  by  a  plate  (see  figure  7)  -which  supported  a  Pirani 
gm  g:,  graasurs  measurement. 

Tb^,  e action  adjacent  to  the  tracking  chamber  was  connected  to  a  mechan¬ 
ical  vacuum  pump  through  a  1  -in.  -diameter  flexible  hose.  Tl.e  two  vacuum 
imca  to  cct.0  differential  pumping  sections  were  supported  by  a  massive  stand 

ia  oVd.r  to  attenuate  vibrations  from  the  mechanical  vacuum  pumps, 

0  0  ^ 

.  t  I 

v  O 

Tha  AGtfc/h'J  .’action  of  the  apparatus  was  the  tracking  chamber.  It  was  .en¬ 
closed  on  Oila  adds  by  the  brass  plate  bearing  the  tracking  system  and  on  the 
sppcsits  cade  by  a  plate  on  which  another  Pirani  gauge,  a  thermometer.,  and 
a  p.r.tccure  £&p  were  mounted  (see  figure  7).  The  bulb  of  the  thermometer 
projected  into  the  chamber  to  sense  the  temperature  of  the  gases. 

o 

"Tb  vpxejcure  tap  was  connected  to  the  manometer  (shown  in  figure  8)  jn 
which  tracldng- chamber  pressure  was  measured.  The  two  working  flu  Is 

in  tb  ^macicmster  legs  were  mercury  and  diffusion  pump  oil,  the  latter  .en¬ 
abling  jSioro  precise  measurements  at  lower  pressures.  The  fluid  levels  were 

'head  'aaihj  the  cathetometer. 


Cac  23  were  fed  into  the  tracking  section  through  the  base  plate.  The 
Saw  of  ip.3  33  and,,  correspondingly,  the  pressure  wer  e  controlled  by  a  fine 
adjustment  of, a  needle  valve  between  the  tracking  chamber  and  the  gas  supply.. 


Cyaicxd  alignment  was  essential  to  the  successful  operation  of  the  expen- 
feoai*  Tha  three  400-micron  orifices  had  to  be  positioned  accurately  on  a 
liras  which  cart-ended  to  the  point  j£  the  magnetized  needle  and  coincided  with 
the  card: crlina  of  the  Faraday  caga  system.  Each  400 -micron  orifice  was 
metre ta-ff  in 'a  larger  plugt which  fitted  loosely  into  3/8-in,  holes  between  the 
£^2.3  sections.  This  loose  fit  permitted  positioning  of  the  orifice  by  lateral 
screws.  Sealing  was  accomplished  through  an  O-ring  face  seal. 
Tbs  olignzQ sht  of  the  orifices  and  the  needle  point  was  carried  out  with  a  light 
esmre Si  < 


2.3  TIt&GSTRS  SYSTEM 
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A”  carles  of  Faraday  cages  was  the  fundamental  .element  of  the  tracking 
,  aysi.im,  As  the. charged  particle  passed  through  each  cage,  its  charge,  or  a 
fraction  tbbr^of,  was  induced  on  the  cage  and  appeared  as  a  voltage  change  on 
the  grid  of  a  .cathode  follower.  Thus,  as  the  particle  passed  each  cage,  a  s^nal 
esuld  b^  observed,  oscilloscope.  ,  '  ' 
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Xhc  Faraday  cages  for  the  velocity -measuring  and  tracking  chambers 
,  were  mounted  ona  single  plate  to  facilitate  alignment  with  the  system,  as 
ehovm  in  figure  9.  The  series  of  cages  on  the  left  fitted  into  the  velocity - 

staring  chamber  while  those  to  the  right  comprised  the  tracking  section. 

o  -  '  < 

Tho  cages  for  the  velocity -measuring  section  consisted  of  a  wire  stub, 
fhree  rings,  and  a  tube.  The  stub,  one  ring,  and  the  tube  are  mounted  on 
a. wire  connected  directly  to  the  grid  of  a  miniature  high-ptriode  (6CW4). 

The  remaining,  avo  rings  were  grounded.  The  signal  generated  by  the  wire 
,  Stub  oa  passage  of  a  particle  served  totrigger  the  oscilloscope. 

0  Q  '  ^  '  c.  * 

Jplevea  rings  made  up  the  tracking- section  cages.  Five  of  these  were 
connected  in  commou  to  the  grid  of  another  6CW4,  and  the  other  six  were 
grounded.  The  cages  were  separated  by  1  -cm  distance  in  ‘he  tracking  section 
.  and  fey  1/2-cm  in  the  velocity-measuring  section.  The  distance  was  measured 
to  within  1/10  of  \%  by  attaching  the  system  to  a  machine  bed  with  a  fine 
perotv  feed  and  noting  when  each  ring,  viewed  by  a  stationary  microscope, 
''reached  the  center  of  the  view  field, 

0.  •»  c 

A  typical  oscilloscope  trace  from  the  Faraday  cage  system  is  shown  in 
figure  iO.  The  scope  was  triggered  by  the  small  wiie  stub;  then  the  signal 
went  through  a  minimum  as  the  particle  passed  the  grounded  cage.  A  maxi¬ 
mum  and  a  minimum  signal  were  produced  as  the  particle  traversed  the  next 
signal  and  grounded  cages,,  and  then  a  longer  duration  signal  was  produced  as 
the  particle  traveled, through  the  tube.  The  length- diameter  ratio  of  this  tube 
pas  so  Targe  that  the  particle  fs  entire  charge  was  induced  on  the  tube;  there¬ 
fore,  the  signal  magnitude  was  directly  proportional  to  the  charge  on  the 
>  v^article..  As  the  particle  proceeded  through  the  tracking  section,  alternating 
'signal  minima  and  maxima  were  produced  as  the  particle  passed  through 
grounded  and  signal  oag.es.f  respectively.  The  increasing  time  interval 
0  between  signal  peaks  corresponds  to  the  particle  deceleration  produced  by 
aerodynamic  drag. '  J- 

b  *C  '  f  - -V 
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The  question  arises  as  to  whether  or  not  the  transfer  of  charge  in  the 

Faradgy  cage  system  may  constitute  a  decelerative  force  on  the  particle. 

This  can.be  evaluated  by  comparing  the  electrical  work  dissipated  to  the 
,totel  change  in  the  kinetic  energy  of  the  particle,  which  is  expressed  by  the 
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where  i  is  current,  R  is  the  resistance  of  the  Faraday  cage  system,  q,hd 
t  is  the  time  for  the  particle  to  pass  through  the  range.  The  current  is  esti¬ 
mated  by 


qu 


where  dc  is  the  distance  between  each  cage.  Using  this  approximation  and 
realizing  that  u^/uq  ~  i,  the  above  ratio  becomes 
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where  is  the  total  range  distance.  The  values  of  the  parameters  are 


q/m  <  1  coulombs /kg 
m  <  4  x  IQ"*1  kg 
d'  =  10-1 
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which  yields 
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R  <  1  ohm 


u  >  10  m/ sec 
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and  demonstrates  that  the  effect  of 'the  electrical  losses  on  the  deceleration 
of  the  particle  is  insignificant. 

As  discussed  above,  a  high-p  miniature  triode  is  employed  in  the  cathode- 
follower  circuits.  The  customary  grid  resistor  was  omitted,  thus  allowing 
a  maximum  input  impedance  ( ~  l6Hf2)  and  power  gain  (~10^  V).  A  low  pass 
filter  is  inserted  in  the  output  to  reduce  false  triggering  of  the  oscilloscope 
by  fast  transients  originating  in  other  equipment.  The  voltage  gain  of  the 
cathode -follower  in  the  velocity-measuring  section  was  measured  to  be  0.953 
to  within  1/2,  of  1%. 

To  interpret  the  voltage  signal  from  the  tubular  cage  in  terms  of  the 
charge  oh  the  particle,  it  is  necessary  to  know  the  capacitance  of  the  system, 
A  special  capacitance  bridge  was  designed  and  developed  to  measure 
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capacitances  as  low  as  4.38  sc  10“2  pF.  The  capacitance  of  the  velocity - 
section  Faraday  eagre  system,  positioned  in  place  and  under  operating  con¬ 
ditions,  was  found  to  be  6.Q0  pF  to  within  1/2  of  1%.  Thus,  the  charge  on  the 
particle  (q)  is  related  to  the  voltage  signal  from  the  tubular  cage  (E^)  by 

s  '  »{  " 

<,  1  Q  =  5.72  E  sc  10  ^  il%  (coulombs)  (8) 
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Two  reference  signals  are  put  on  each  oscilloscope  trace  for  measure- 
meat  pa&pcrjfce.  A  signal  from  the  internal  voltage  calibration  of  the  oscillo- 
ccopo,  v.'juicix  Xiao  calibrated  against  a  standard  cell,  provides  a  reference  for 
measurement  of  the  Signal  magnitude  from  the -tubular  cage.  A  trace  from-  a 
time-mark  generator  (Tektronix  180A)  serves  to  reference  the  instant  of 
particle  passage  through  each  cage.  These  calibration  signals  are  also  shown 
on  figure  10. 


3.0  OPERATION  OF  THE  APPARATUS 


.  v» 

Ihe  aim  of  the  experiment  was  *6  obtain  drag  coefficient  data  over  the 
entire  M-Re  regime  encountered  by  particles  in  a  rocket  nozzle.  The  vidua  of 
M  is  determined  by  the  particle  velocity  and  the  speed  of  sound  in  the  tracking  - 
section  gas.  The  value  of  Re  is  a  function  of  the  particle's  velocity  and  sites, 

the  pressure  in  the  tracking  section,  and  the  working  gas, 

‘  *  »  ■‘N 

Each  series  of  experiments  was  preceded  by  a  thorough  cleaning  of  the 
acceleration  chamber  and  replenishing  of  the  particle  supply.  After  re¬ 
assembly  the  tracking  section  was  flushed  and  the  desired  pressure  set  by  a 
careful  adjustment  of  the  needle  valve.  A  period  of  approximately  30  min 
was  required  to  stabilize  the  pressure.  The  electronic  instrumentation  was 
then  set  and  the  accelerating  potential  applied  to  the  system.  The  micrometer 
needle  was  screwed  down  by  a  long,  plastic  red  attached  to  the.  micrometer 
handle  until  the  most,  favorable  location  was  reached  and  particles  of  the  de¬ 
sired  charge  and  velocity  were  observed. 

The  signals  from  the  cathode-followers  were  fed  into  a  Tektronix  549 
storage  oscilloscope,  and  undesirable  signals  were  erased.  When  a  particle 
of  the  desired  velocity  and  charge  was  observed,  the  oscilloscope  intensity 
was  quickly  turned  down  to  prevent  additional  traces  from  appearing  on  the 
screen.  An  electronic  circuit  was  also  developed  to  translate  the  base  line 
of  each  new  signal  a  small  amount  vertically  so  that  the  base  lines  of  several 
traces  could  be  distinguished. 

»  '  >'*  'c 

After  selecting  and  isolating  the  desired  signal,  the  accelerating,  voltage, 
tracking-chamber  pressure,  and  temperature  were  recorded.  Then  the  volt¬ 
age  and  time  calibrations  were  put  on  the  oscilloscope  screen.  An  X-Y  plotter, 
adapted  to  the  oscilloscope,  was  used  to  write  an  identifying  number  on  the 
screen.  Polaroid  and  35 -mm  photographs  were  then  taken  of  the  trace. 

Operation  of  the  experiment  is  illustrated  in  figure  li„  The  operator  is 
turning  the  micrometer  drive  to  feed  more  particles  into  the  accelerating 
chamber  while  prepared,  with  his  left  hand,  to  erase  undesired  traces. 

Usually  a  series  of  experiments  were  conducted  at  one  pressure  level  by 
selecting  velocities  ranging  from  the  smallest  to  the  largest  possible.  The 
pressure  range  was  established  by  the  observed  velocity  decay.  If  the  pres¬ 
sure  was  too  low,  insufficient  decay,  occurred  for  accurate  deceleration 
determination;  and  if  the  pressure  was  too  high,  the  velocity  change  became 
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too  large  to  represent  the  test  by  single  values  of  M  and  Re.  The  velocity 
decays  for  the  data  reported  here  varied  between  10%  and  40%.  After  each '  ‘  ' 
test  the  approximate  M  and  Re  tve.  --s  plotted  to  note  the  relation  of 
numbe*.  to  the  flow  regime  of  ?n  esreat. 

Nitrogen  and  Freon  114  were  employed  as  working  gases.  The  h x<y:>2:c‘ 
mo’eculai  weight  of  Freon  (171  g/rifOle)  corresponds  to  a  lower  speed  >;  C 
sound  <~nd  enabled  highe**  M  to  be  achieved.  The  pressure  in  the  tracing 
section  was  varied  between  1  and  50  terr  when  nitr'gen  was  used  and  batWiiva 
1  and  10  torr  with  Freon.  .  e 

°  o 

Two  kinds  of  magnetic  particles,  carbonyl  iron  and  nickel,  were„w s«&d  in 
the  experiments,  A  scanning  electron  micrograph  of  the  carbonyl-iron 
particles  is  shown  ir.  figure  12.  It  is  evident  that  the  particles  were  not 
perfect  spheres;  t.ome  appear  to  have  blisters.  However  the  data  obtained 
should  still  be  indicative  of  the  drag  coefficients  of  particles  in  a  rocket 
nozzle.  The  nickel  particles,  on  the  other  hand,  were  quite  spherical  as 
illustrated  in  figure  13,  The  size  of  carbonyl-iron  particles  measured  in 
the  tests  ranged  from  1. 6  microns  to  6  microns  and  the  nickel  particles 
from  5  microns  to  17  microns  in  diameter.  This  agrees  reasonably  with 
the  sizes  observed  in  the  photomicrograph.  \  g  ' 

o'  *_  '  a 

The  range  of  M-Re  covered. >-y  the  experiment  is  illustrated  in  figure  14. 
A  flow  regime  which  extended  from  free-molecule  flow  through  transition  * 
and  just  into  slip  flow  was  investigated.  Thus,  the  trajectory  of  a  particle 
in  a  typical  rdeket  nozzle  lies  within  the  range  covered. 
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Figure  13.  Electron-Microscope  Photograph 

of  Nickel  Bar  tides 
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4.  0  DATA  REDUCTION 


The  fundamental  parameters  to  be  reduced  from  the  data  were  tha  par iic:L.z 
drag  coefficient  and  the  corresponding  id  and  Re.  The  following  sactio&s  dip-  e 
cuss  the  data  reduction  technique,  its  accuracy,  and  its  application  to  a  typical 
test. 

4. 1  CALCULATION  OF  THE  DRAG  COEFFICIENT  FROM  DATA 

£>  <  O  < 

The  aerodynamic  drag  force  on  a  particle  is  equal  to  the  product  of  its  « 
mass  and  deceleration.  The  particle  mass  is  determined  when  it  passes 
through  the  velocity- measuring  section,  and  its  deceleration  is  determined 
from  the  velocity  decay  in  the  tracking  section. 

C  c  ' 

Twp  differentiations  of  time-  distance  data  to  deduce  deceleration  can  w 
lead  to  significant  errors,  it  is  more  desirable  to  use  an  expression  which 
directly  relates  drag  coefficient  and  the  time- distance  data.  Equating  the 
aerodynamic  drag  on  a  particle  and  its  deceleration  yields* 
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Over  the  range  of  velocity  decay  in  these  experiments,  it  was  reasonable  to 
assume  that  the  drag  coefficient  was,  in  essence,  constant.  The  above  equa¬ 
tion  then  can  be  integrated  to  give  the  following  relation  fox  the  time- distance 
history  of  the  particle  - 
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where  uD  is  the  initial  velocity  of  the  partible,  s  the  distance  traveled  in 
time  t,  and  t0  the  time  when  s  =  o  .  a  is  a  constant  related  to  the  drag 
coefficient  by 
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*It  can  be  shown  that  the  acceleration  due  to  gravity  can  he  neglected  for  the 
range  of  particle  decelerations  of  these  tests. 
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where  Pp  and  rp  are  the  particle  density  and  radius,  respectively.  It  is 
convenient  to  introduce  a  new  variable  t)  where 
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6  Amplifies  to 
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For  a  given  a  ,  the  least- squares  technique  can  be  used  to  determine  values 
of  and  to,  corresponding  to  a  straight  line  or.  a  plot  of  t  versus  rj  , 
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is  a  minimum.  oThis  was,  repeated  for  deferent  values  of  a  to  find  tha.  value 
c-l  c  for  w^ich  E(or)  is  a  minimum .  (Because  the  spatial  location  of  each 
cage  was  measured  to  0, 1%  there  was  no  need  to  consider  deviations  in  s  . ) 
Having  determined  the  va,l,ue  of  a  which  represented  the  best  fit  of  the  data, 

fhc' drag  coefficient  was  determined  from 
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A  computer  program  was  developed  to  calculate  for  the  minimum  squared 
deviation  of  the  data.  *  c 
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An  experimental  determination  of  the  particle  density,  using  standard 
aspdytical  techniques,  showed  that 
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o  °  (carbonyl  iron)  =  7.26  g/cm 


p  (nickel);,  c<  =  8.52  g/cm" 
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to  "within  1/10  of  1%.  <•-'  v 
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The  gas  density  was  calculated  from  the  temperature  and  pressure 
msacuromeats  in  the  tracking  chamber,  using  the  equation  of  state  for  an 
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The  remaining  parameter,  rp ,  was  found  from  the  data  obtained  by 
means  of  the  velocity- measuring  section.  Representing  the  particles  a3 
spheres,  one  can  rewrite  equation  3  as 
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The  measurement  of  q,  V,  and  the  particle  velocity,  u_f,  has  been  dis¬ 
cussed  above. 

Values  of  M  and  Re  were  calculated  using  the  average  particle  velocity, 
i.e.,  the  distance  between  the  first  and  last  tracking  cages  divided  by  the 
transit  time. 

4.  2  ERROR  ANALYSIS 

The  accuracy  of  the  calculated  drag  coefficient  depends  on  the  accuracy 
with  which  Pp ,  Pg ,  rp=,  and  a.  can  be  determined.  The  variance  (cr4)  of 
the  drag  coefficient  is  related  to  that  of  the  independent  variables  through 
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where  the  contribution  due  to  p^  has  been  neglected  because  of  the  high  ac¬ 
curacy  with  which  it  was  determined.  Following  standard  statistical  proce¬ 
dures,  the  variance,  or  standard  daviation,  can  be  related  to  the  magnitude 
of  error  at  a  chosen  confidence  level. 

The  variance  of  a  is  determinable  from  the  tracking- section  data.  It 
is  related  to  the  variance  of  the  time  deviation  ^<7^  j  about  the  least  square 


curve  by 


a 


■  £  W- 


(18) 


K 


N> 

O 


i  O0 


o 

t 

o* 

| 

4 

4 


^  "*> 


A  seaming  the  variance  is  the  same  for  all  points  yields 
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:czn  be  ohowrx  that  the  sum  of  the  derivatives  of  a  with  respect  to  t^  , 
cebjcct  to  the  constraint  of  a  minimum  sum  of  squared  deviations  about  the 
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the  subscripts  K  refer  to  the  various  distance  stations,  and  N  is  the  number 
of  stations  in  the  tracking  section. 
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The  best  estimate  of  the  variance  of  the  time  deviations  is 
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where  N  is  the  number  of  stations  in  the  velocity- measuring  section.  For 

ibe  of  the  present  tests  it  was  found  that 
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sal  tho -error  is  r  could  be  neglected, 

<•  5  •  P  .  -  •  •  -  o 

Th®;  variance  of  the  gas  density  is  related  to  the  uncertainty  in  pressure 

ated  tos^porature  ssiea sur  ©meats  by 


,Oseo  again  it  was  found  that  this  variance  is  much  smaller  than  f—2L  )  and 

*  o  <  ..  ,  c  ,  '•  '  a  ' 

<could5bs  neglected, ,  -  «  v  ?  ' 
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-  //  Thus,  the  variance  of  the  drag  coefficient  becomes  simply 
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and5 tlie  error  is- 
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where  "Q  wae-determinedfrom  the  Student  -t  iistribution  for  the  75%  con- 
iidettce  level, 

,  i.  *  C 

'  *0  o  - 

,  0  w* 

4.5  .  REDUCTION  OF  A  DATA  POINT  FROM  A  TYPICAL  TEST 

.  ,  ,  *  ,  c.  ' 

The  trace  from  a  typical  test,  using  carbonyl-iron  particles  in  nitro¬ 
gen,  is  shown  in  figure  15.  The  small  time  divisions  represent  50  psec, 
the  larger  divisions  0.  5  msec  intervals.  The' volt  age  calibration  repre¬ 
sents  10  mv*  The  pressure  for  this  test  was  measured  as  2.  24  cm  of 
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dlfjfecica  p-nmp  ©H>  the  temperature  was  2S.5°C,  and  the  potential  on  the 

26, 0  volts. 
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Tha  f.kspt  parameter  calculated  was  the  particle  radius  which  was  deter- 

•  nJ'  \  IJ.sk  -A  the  data  o t  the  velocity  measuring  section.  In  measuring  the  die* 

^ trace  between  the  two  grounded  cages  and  relating  it  to  the  time 

co*-  ;<v  y.c-ra,  £t  to  fennd  that  the  particle  traveled  1  cm  in  70  peoc  or  had  a 

of  10  m/a@c.  Multiplying  the  electrometer  voltage  by  the  calibration 

iacln^  for  the  voltage  divider,  1,9.81  x  10  ,  gives  an  acceleration  potential 

of  SI, Si  ?.:v.  The  magnitude  of  the  potential  induced  on  the  tubular  Faraday 

cage  io  measured  as  li,35  mv  from  the  trace  and  the  particle  charge  calcu- 

“  '  - 14 

lak  $  from  equation  4  is  q  =  6,5  x  10  coulombs.  Substitution  of  these 

data.  and  the,  donaity  ©f  carbonyl  iron  into  equation  16  yields 

t  <  o  '  0 

+  ,  o  '  * 

'  4  ‘  r  =  2,2  microns 

,  i.  %  -■  "  p 

o  ;  ^ 
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The  gao,  density  is  derived  from  the  pressure  and  temperature  measure- 
moats,  The  density  of  the  diffusion  pump  oil  was  measured  to  be  0. 9825  g/ 

,  a  ,  ■  ”  3 

era  and,  therefore,  2,  24  cm  of  oil  is  equivalent  to  1.62  torr,  or  2,  16  x  10 

2 '  -  '  ‘ 

dyne  s /cm  .  This  value  and  a  temperature  of  25.5°C  give 

© 

,  '  P  2.44  x  10~  g/cm 

„  «  '  g "  ,  ° 

-  s  -v*  '■  -  °  '  ' 

Feeding  the  time-distance  data  measured  from  the  tracking- section 
signal  trace  into  the  computer  program  to  determine  a,  and  the  correspond¬ 
ing  leads  to  ..  . 


,  *  -  ,  a  =  0.00858  *11% 
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at' a  ?S%  confidence  level.0* 
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,  »  Substituting  the  above  values  into  equation  15  results  in 


'  O  o  ^  o 


CD  =  15.  Oil  1% 


The  .average  velocity  in  the  tracking  section  was  found  to  have  boon 
134  m/ sec  which  yields  I  1  „  -  -  . 

,,  c  i 

<e  ,  /  P.e  =  0.081  and  M  =  0.38 
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The  corresponding  Kn  is  7. 1,  and  is  related  to  M  and  Re  by 

Kn  =  1.  26  \[y  M/Re 

After  completion  of  the  experiments,  it  was  discovered  that  a  leak  had 
developed  in  the  tracking  section  and  the  lower -pressure  Freon  results  were 
appreciably  affected.  The  magnitude  of  the  leak  was  determined  by  comparing 
nitrogen  and  F. >pn  data,  which  overlapped  in  certain  flow  regimes,  and  all  the 
Freon  data  were  corrected  accordingly.  c 

A  summary  of  the  data  is  given  in  table  I. 

TABLE  I 

SUMMARY  OF  EXPERIMENTAL  DATA 
Particle  Type  Particle 


Test  No. 

and  Gas 

M 

Re 

Kn 

Drag  Coefficient 

1 

Carbonyl- Iron 

0,65 

0.36 

2.7 

7.4*12% 

2 

Nitrogen 

0.88 

0.27 

4.8 

5.  5*15% 

3 

0.97 

0.18 

8.0 

6.3*16% 

4 

0.25 

0.09 

4.1 

17*13% 

5 

0. 18 

0.042 

6.4 

28*6% 

6 

0. 11 

0.032 

5.1 

52*6% 

7 

0.061 

0.014 

6.5 

90*22% 

8 

0.45 

0.19 

3.5 

9.9*14% 

9 

0.21 

0.22 

1.4  ■ 

2-1*11% 

10 

0.37 

0.45 

1.2 

11*9% 

11 

0.  16 

0. 14 

1.7 

27*14% 

12 

0.  11 

0.010 

16 

61*11% 

13 

0. 18 

0.067 

3.9 

30*10% 

14 

0.62 

0. 11 

8.4 

8. 3*7* 

15 

0.53 

0. 11 

7.2 

12*13% 

16 

0.84 

0. 13 

9.6 

7.7*17% 

17 

0.38 

0.080 

7.1 

15*11% 

18 

0.30 

0.045 

9.9 

18*16% 

19 

0.48 

0.15 

4.8 

9.9*11% 

20 

0. 13 

0.071 

2.7 

38*13% 

21 

0.25 

Q.  13 

2.8 

21*10% 

22 

0.82 

0.21 

5.8 

6.8*2% 

23 

0.092 

0.075 

1.8 

62*12% 

24 

Carbonyl- Iron 

0.060 

0.032 

2.8 

79*9% 

25 

Nitrogen 

0.28 

0,020 

21.4 

21*8% 
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SUMMARY  OF  EXPERIMENTAL  DATA  (Continued) 

*  1 

•_ '  Particle  Type 


Particle 


■’  ‘  v^-' '  \ V,  * 

*  «?v  «.  fl,  >.e  '  % 

aa4  Gas 

M 

Re 

Kn 

Drag  Coe£fici< 

0  ** 

\  *  e  ■ 

"  “ 

-  *s 

Carbonyl -Iron, 

-  0.65 

0.33 

2.6 

8.0*12% 

t*  #  0  ; 

Freon 

1.48 

0.97 

2.0 

3.7*12% 

x> 

22  .  « 

c 

v  y 

1.77 

0.69 

3.4 

4.0*11% 

«a>  ‘  s  < 

29  " 

"<  4 

o 

0.97 

0.23 

5.5 

5.8*8% 

30, 

0.78 

0.41 

2.5 

7.8±127o 

'■  ?1  " 

0.96 

0.82 

1.55 

5.3*7% 

>> 

-.32.  / 

<  Carbonyl-Iron 

1.3 

1.20 

1.42 

4. 9*4% 

ft 

=-  33 

..  Freon 

1.22 
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5.0  DATA  ANALYSIS  .  , 

c 

(7)  (8  9  10) 

Sherman'  '  has  shown  that  drag  coefficient  data'  *  ’  -of  sphere,?,  c.t 

high  M  (M  >  4)  can  be  reduced  to  a  nondimensional  drag  coefficient  which  , 

depends  only  on  Kn,  i.  e. ,  the  ratic  of  the  mean-free  path  of  the  gas  moloc#tor> 

to  the  particle  diameter.  The  nondimensional  coefficient  i3 

C  -  C 
D  D, 


where  C  *  is  the  drag  coefficient  for  Re  -*  <»  and  C 

°FM 

molecule- flow  value. 


the  free- 


Ballistic-  range  data'  '  at  high  Re  indicate  a  constant  drag  coefficient  of 
0.  92  for  M  greater  than  2.  This  value  was  used  jv,  •  C  in  this  M  range. 

DI 

However  no  data  are  available  for  the  drag  coefficient  at  extremely  large  Re 

7 

(Re  >  10  )  and  low  M.  For  purposes  of  the  present  atv  lysis,  the  variation 

5 

of  drag  coefficient  with  M  at  Re  ~  10  ,  shown  in  figure..  *6*  will  be  used 
for  the  C  in  equation  30. 

T-he  drag  coefficient  of  a  sphere  in  free-molecule  flow,  assuming  diffuse 
reflection  of  the  molecules,  is  given  by  ^  ‘ 

c  =  (1  +  2S2)  +  Is— t-i|- 1  erf  (S)  +  4~-(31) 

FM  S  2S  5  '  % 

where  S  is  the  speed  ratio,  which  is  related  to  M  by 


*Sherman  refers  to  this  parameter  as  the  "inviscid"  drag  coefficient. 
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The  above  equation  assumes  thermal  equilibrium  between  the  particles  cx*4.  ^ 
the  gas.  For  low  M  it  can  be  shown  that  .  ■. 


FM 


6,3 

•JyU 


•(M  <<  1) 


*$4 


Drag  coefficient  data  obtained  in  the  present  experimental  studies  h'vv-a. 
been  reduced  to  the  nondimensional  parameter  given  by  equation  30  an.d  are  ,,,  v. 
plotted  versus  Kn  in  figure  17,  Data  from  other  experiments  at  high  v 

^8  9  jo)  (12}  '  .c * '' 

M,  ’  ’  '  Millikan's  results'  1  from  the  classical  oil  drop  experiment,  .and 

(13)  ^ 

low  Re  data  collected  with  a  magnetic  suspension  system  are  shown  on  the 
same  graph.  At  high  Kn,  the  data  converge  to  one  value,  the  free  molecule-  • 
flow  limit.  Also,  for  decreasing  Kn  the  data  for  high  M  tend  to  group  around 
a  single  curve.  However,  the  drag  coefficients  corresponding  to  low  Re  tend 
to  break  away  from  the  high  M  limit  with  decreasing  Kn.  , 

At  very  low  Re  in  continuum  flow  where  Stokes*  law  is  valid,  the 
drag  coefficient  is 


4 


24 

*mmum** 

Re 


(34)  c 


<r  c 
by  O 


Correspondingly,  the  free- molecule- flow  result  is  given  by  equation  33  and 
the  nondimensional  drag  coefficient  approaches 

S*  o^Sr151315”  <3S> 

as  Kn  -*•  0 ,  Thus,  at  low  Re ,  the  curve  must  break  away  to  become 
asymptotic  to  the  above  relation. 


Similarly,  the  necessary  asymtote  at  low  Kn  for  curves  at  higher  Re  is 

°*  / 

(36) 


Gj^  =  Kn  Re 


(s,  -  °-48) 

\  me _ / 


C\- 


8 


where  is  the  drag  coefficient  for  incompressible  flow, 

inc 

No  theoretical  analysis  is  available  which  provides  an  analytic  expression 
for  the  variation  of  with  Kn  and  Re  .  Some  analyses ^  have  been 

published  assuming  near- free  molecule  flow,  but  these  are  valid  only  for  Kn 
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of  unity  and  greater.  A  complex  Monte-Carlo  technique  has  been  developed 

o  /  i  4\ 

by  Vo£eni$J5,  et  al.,  '  to  treat  flows  from  continuum  to  free- molecule  flow. 
This  technique  cuececsfully  predicts  the  measured  drag  coefficients  at  high 
"M  but  docs  not  yield  an  analytic  expression  for  versus  Kn, 

v  ~  K  ,  " 

'*  It  is  interesting  to  note  that  Vogenitz,  et.  al.,  predict  a  drag  coefficient 
at-  Kn  =  100  somevyhat  higher  than  the  analytic  result  calculated  from  the 
conventional  free- molecule- flow  analysis.  The  same  tendency  is  observed  in 
th®  results  of  the  present  study.  Millikan’s  results  fall  somewhat  below  the 
predicted  Cn  ,  which  is  probably  attributable  to  specular  reflection  from 
<  FM  * 

the.  smooth  oil  drops .  ?  *  - 

°o: 

Axi  empirical  expression  which  reasonably  fits  the  high  M  data  .s 


M»1  t 


1  1.1  Kn  ^*3  exp  (-Kn*^) 


=  G(Kn) 


and  is  shown  in  figure  17.  It  is  assumed  that  the  C  approaches  1. 1  as  the 

<  Kn  extends  to  infinity.  r*  ' 

/  '  .  *'•  *  *  '  •  ° 

%  c  Q  °  ,  ' 

o  In  order  to  obtain  an  equation  representing  all  of  the  data,  equation  37 
can  be  multiplied  by  an  empirical  factor  expressing  the  break- away  of  the 
low  Rc  curves.  In  the  low  Kn  region,  the  function  G(Kn)  is  closely  repre- 

'  <>..  sealed  by  *  ‘  v 


•  "  ,  G(ICn)  UJ  O.SE  Kn0*3  (i0“3  <  Kn  <  10rl)  ( 

v  O  ° 

By  comparing  equations  3*  and  36,  it  is  seen  that  at  low  Kn  the  empirical 
factor  must  approach? 


„o,7(\ 


0 . 48 1  Re 


to  give  the  correct  asymptote.  A  factor  incorporating  this  feature  and  provid¬ 
ing  an  acceptable  fit  with  the  data  is 


i>'Kn,  Re)  =  j  1  -  exp  J-Kn  e 


Kn  \_inc 


0.48^  Re 

>.6 


Figure  17.  Nondimensional  Drag  Co-  ( 
efficient  for  a  Sphere 
vs  .Knuds  en  Number 
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Thus,  the  complete  empirical  equation  for  is 

CD  =  G(Kn)  •  D  (Kn,  Re)  (40) 

Empirical  curves  for  Re  =  200,  10,  and  less  than  unity,  calculated 
from  equation  40,  are  shown  in  figure  17.  Agreement  with  the  available  data 
is  good  and  equation  40  is  recommended  for  calculating  the  drag  coefficients 
of  particles  in  a  rocket  nozzle. 


It  is  enlightening  to  trace  the  trajectory  of  a  particle  in  a  rocket  nozzle 
with  reference  to  figure  17.  The  Kn  of  a  particle  is  given  by 


Kn  =  1.26 


47 \i. 
g 

2r  p  a 
P  g 


where  jx  is  the  gas  viscosity  and  a  the  local  speed  of  sound.  Because  the 
g 

gas  density  decreases  monotonicaily  with  distance  through  the  nozzle,  the  Kn 
monotonically  increases.  The  Kn  of  a  3- micron  particle  in  a  200- psi  rocket 
chamber  is  0.031.  Thus,  the  trajectory  begins  on  the  Re  <  1  line  at  the 
chamber  Kn ,  goes  through  the  maximum  Re  in  the  throat  region  where 
Kn  ~  0.06,  and  proceeds  back  toward  the  Re  <  1  curve  with  increasing 
Kn ,  traversing  the  whole  flow  regime  covered  by  this  study. 

It  is  interesting  to  compare  the  data  with  expressions  devised  for  drag 

(2) 

coefficient  versus  M  and  _le  before  the  data  were  available.  Crowe  '  pro¬ 
posed 

0 

CD  =  (CD.  ”  *).«*{-'».« Jr  £*  0W)1  +  *  Z  (42) 

\  me  V  I  i  7  M 


where  lqg1Q  gi  (Re)  =  1.25  Jl  +  tanh  (0*77  log1Q  Re  -  1.92) 


and  ' 


h(M)  *  2.3  +  1.7  (T  /T  ) 

:P  g 


?  -  2. 


3  tanh  (1.17  log1Q  M) 


OH  °  0  ,  a  y 


S’ 


-  <*C 


C  c  <■  " 


T  /T  is  the  particle  temperature  —  gas  temperature  ratio  which,  for  this 

P  ®  (17) 

comparison,  is  unity.  Carlson  suggested  the  formula 


C 


(1  +  0.15  Re0*687) 

i  1  +  exp(-0.427/M^*^8)  exp  (-3/Re^‘88)] 

‘  M 

1  +  Rll 

3.82  +  1.28  exp  (-1.25  Re/M), 

(43) 


/  1Q\ 

based  on  various  analytic  and  empirical  trends.  Kliegel'  ,  in  his  two- 
phase  flow  analysis,  has  used  the  drag  coefficient  predicted  by  solving  the 
13  moment  equations  for  flow  over  a  sphere,  namely 


C 


D 


xnc 


(i  + 


(1  +  7. 5  Kn)(  1  +  2Kn)  +  1.91  Kn2 
7.  5  Kn)(  1  +  3  Kn)  +  (2.29  +  5. 16  Kn)  Kn2 


(44) 


The  above  three  expressions  for  drag  coefficient  are  plotted  on  figure  18. 
The  various  data  are  also  shown  for  reference.  All  the  curves  correspond  to 
an  Re  of  unity.  It  is  noted  that  all  the  curves  become,  parallel  with  the 
asymptotic  line  for  Stokes  flow  at  low  Kn  (C^  “3  Kn),  but  are  somewhat 

higher  because  the  drag  coefficient  at  Re  =  1  is  slightly  larger  than  pre¬ 
dicted  by  Stokes  drag  law.  With  increasing  Kn,  ithe  drag  coefficient  expres¬ 
sions  used  by  Carlson  and  Kliegel  diverge  significantly  from  the  data*  The 
expression  proposed  by  Crowe  shows  better  agreement  but  displays  an  unde¬ 
sirable  inflection  in  the  near- free- molecule  flow  regime,  near  Kn  =  6  . 

Such  an  inflection  is  not  predicted  by  present  theories,  and  the  use  of 
equation  42  is  not  recommended. 
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Figure  18.  Comparison  of  Proposed  Expressions 
for  Particle  Drag  Coefficient 
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6.0  CONCLUSION 


Drag  coefficients  for  spherical  particles  were  measured  over  a  range  of 
M  and  Re  corresponding  to  flow  regimes  encountered  by  particles  in  a  rocket 
nozzle.  The  data  were  reduced  in  terms  of  a  nondimens ional  drag  coefficient 
and  correlated,  together  with  available  data  obtained  under  other  flow  condi¬ 
tions,  in  terms  of  two  fundamental  parameters:  Re  and  Kn.  An  empirical 
equation,  valid  for  all  Re  and  for  Kn  from  0.  Q01  to  infinity,  was  developed 
for  the  nondimensional  drag  coefficient.  This  equation  is  recommended  for 
use  in  calculations  of  gas-particle  flow  in  rocket  nozzles. 
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